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Abstract 

We discuss here in detail a new analytical random walk approach to calculating the 
phase-diagram for spatially extended systems with multiplicative noise. We use the 
Anderson localization problem as an example. The transition from delocalized to 
localized states is treated as a generalized diffusion with a noise-induced first-order 
phase transition. The generalized diffusion manifests itself in the divergence of aver- 
ages of wavefunctions (correlators). This divergence is controlled by the Lyapunov 
exponent 7, which is the inverse of the localization length, ^ = I/7. The appear- 
ance of the generalized diffusion arises due to the instability of a fundamental mode 
corresponding to correlators. The generalized diffusion can be described in terms of 
signal theory, which operates with the concepts of input and output signals and the 
filter function. Delocalized states correspond to bounded output signals, and local- 
ized states to unbounded output signals, respectively. The transition from bounded 
to unbounded signals is defined uniquely by the filter function II{z). 
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1 Introduction 



A variety of phenomena in pliysics, chemistry and biology is modelled by 
stochastic differential (or difference) equations [1,2] which involve both addi- 
tive and multiplicative noise. The additive noise is relatively easy to handle 
with the help of the central limit theorem. This type of noise disappears from 
the noise-averaged form of the equations. The situation changes dramatically 
with the appearance of multiplicative noise [3]. Some situations have been 
documented in the literature, in which the multiplicative noise participates in 
the creation of a pure noise-induced phase transition [4]. In such funda- 
mental theorems of an importance comparable to the central limit theorem are 
still lacking. In this paper we will investigate another famous example of this 
phenomenon - the Anderson localization as a noise-induced metal-insulator 
phase transition. 

Let us consider spatially extended systems with multiplicative noise. Anderson 
localization [5] has remained a hot topic in the physics of disordered systems 
for a long time (see review articles [6,7,8]) and constitutes a multidisciplinary 
problem. Formally the Schrodinger equation with random on-site potentials in 
the tight-binding representation is a stochastic algebraic equation with multi- 
plicative noise, where physical meaning can only be attributed to certain aver- 
age values (averaging over the ensemble of random potentials [9]). The prob- 
lem is very similar to statistical physics, especially to the statistical physics 
of phase transitions (quantum phase transitions, see [10]), because there also 
a noise- induced metal-insulator phase transition is analyzed (keywords: crit- 
ical dimension, correlation length, etc.). In the case of phase transitions one 
finds that the relevant parameters (in our case e.g. Lyapunov exponents or lo- 
calization lengths) are not analytical functions of the disorder. Typically the 
theoretical analysis of disordered systems is based either on approximations or 
on numerical methods [6,7]. However, it is clear that this cannot be sufficient 
for systems revealing a phase transition [11], e.g. a metal-insulator transition 
as is the case with Anderson localization. In such a case an exact solution is 
greatly needed. 

Disordered systems are usually considered as a part of solid state physics (or 
statistical physics). The theoretical methods here are very different from the 
ones in the theory of dynamical systems. It is known that static randomness 
(multiplicative noise in space) in the physics of disordered systems may dras- 
tically change macroscopic quantities. Such an effect exists also for dynamical 
systems with randomness in the equations of motion (multiplicative noise in 
time). For example, in special cases the random walk problem is mapped to 
a quantum mechanical tight-binding model exhibiting Anderson localization 
[12,13,14]. It might be reasonable to start from a simple microscopic stochastic 
kinematic model (Brownian motion and normal diffusion), and to build the 
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theory of Anderson localization along this line (tight-binding model and gen- 
eralized diffusion [13], with a noise- induced metal-insulator phase transition). 

In our papers [12,13] we presented somewhat schematically and by leaving out 
many mathematical details an exact analytic solution of the Anderson local- 
ization problem. Subsequently there were critical comments [15] (see also [16]) 
that our proof is too short and condensed, and more details arc needed, the 
more so, since the mathematical formalism used by us is quite new for this 
scientific community of Anderson localization. This is why we present in the 
present paper the detailed derivation of the analytical solution with illustra- 
tions. Note that under the exact solution we mean the calculation of the phase 
diagram for the metal-insulator system [12,13]. We do not calculate transport 
and other important properties. However, the knowledge of the phase diagram 
permits to understand, how these properties can be calculated. 

The structure of the present paper is as follows. In Section 2 we explain how 
the localized states can be treated in terms of a generalized diffusion. This 
approach allows to understand why for defining the phase diagram it is suffi- 
cient to solve exactly only equations for the joint correlators. We present and 
solve here the equations for arbitrary dimensions. Section 3 deals with the 
stability of this solution. It is shown that the generalized diffusion arises due 
to the divergence of the fundamental mode. The determination of the stability 
range of this mode permits us to calculate the phase diagram. We explain 
our language of input and output signals, filter function which is new for the 
Anderson localization community. 



2 Equations for correlators 

2.1 General aspects 

Strictly speaking, the disordered Anderson tight-binding model [5] with Schro- 
dinger equation 

^^M' = {E - eM)^M, (1) 

M' 

where the summation over A4' runs over the nearest neighbours of site — 
{mi, 7712, ■ ■ ■ jITT'd}, cannot be solved exactly for arbitrary dimension D. The 

reason is that the random potential £m enters the equation as a product with 
the random amplitude ipM which corresponds to the multiplicative noise case 
[12]. Therefore, the exact solution of eq. (1) discussed in [12,13] is possible 
only under special circumstances to be discussed below. 
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The fundamental quantity of a disordered system - the locahzation length ^ - 
was determined in [12,13] via the Lyapunov exponent 7. In these calculations 
the following conditions which appear contradictory at first sight have to be 
satisfied. The phase diagram of the system with metal-insulator transition 
should be obtained in the thermodynamical limit [11] (the infinite system). 
However, the determination of the Lyapunov exponent needs the introduction 
of a coordinate system (starting point) which imposes certain limitations on 
the system size. Moreover, a direction for the growth of the divergent quantity 
should be chosen, despite the fact that all space directions are equivalent in 
eq. (1). 

All these conditions are fulfilled for the semi-infinite system, or system with 
a boundary, where the index n = m,D > 0, but all m,j G (— cxd,cxd), j ~ 
1, 2, . . . ,p, with p — D — 1. The boundary which is the layer n = defines the 
preferential direction (the axis n) along which the Lyapunov exponent 7 will 
be calculated. 

It is convenient to interpret the index n — 0, 1, . . . , 00 as the discrete-time, 
whereas all other indices combine in the vector m = {^i, 777.2, ■ ■ ■ , "7p}. The 
Schrodinger equation (1) can be rewritten as a recursion equation (in terms 
of the discrete-time, and assuming summation over repeated indices) 

'0n+l,m £n,m'0n,m V'n— l,m ~l~ ^m,m'^n,m' i (2) 



where the operator 



(3) 



is introduced for compactness (summation over m" includes the nearest neigh- 
bours of the site m) . The discrete-time equation (2) is a difference equation 
of the second order, which needs two initial conditions (Cauchy problem for 
the difference equation). The first natural condition is ipo,m = 0. The second 
initial condition can be presented in the general form ipi^m = without 
any additional conditions for the arbitrary constants am, except that they 
are supposed to be finite. The rules for the treatment of the field am will be 
explained below, section 3.3. 

The recursion equation (2) reveals a general feature of the causality - in its 
formal solution ipn+i,m depends only on the random variables £„',m' with n' < 
n. This gives us a hint for the exact solution of the equation provided that 
random variables Em on different sites do not correlate: (em^M') oc 6m,M'- In 
this case all amplitudes ip on the r.h.s. of eq. (2) are statistically independent 
of £n,m- In other words, while performing the mathematical operations in eq. 
(2), e.g. taking the square of both equation sides, with a further averaging 
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over an ensemble of different realizations of the random potentials (symbol 
(...)), the average of the product of amplitudes iJj and potentials 

^n.m can 

be replaced by a product of the corresponding average quantities. Therefore, 
the exact solution can be obtained only when the on-site potentials Sn^m are 
independently and identically distributed. We assume hereafter existence of 
the two first moments, (en,m) ~ and (^^^m) — "^^^ where the parameter a 
characterizes the disorder level. 

2.2 Generalized diffusion 

The ensemble averaging using recursion eq.(2) allows us to calculate averages 
of different functions containing the amplitudes ip. The problem arises, which 
averages we can and should calculate? E.g. several even moments of the ampli- 
tude were calculated in the 1-D case [17]. whereas only the second momentum 
was calculated for multidimensional systems in [12,13]. It is generally believed 
that the complete information is available from the complete set of all mo- 
ments [18]. However, it is not clear, how this information should be analyzed. 
In particular, how can the one parameter of our interest - the localization 
length ^ - be extracted from the infinite set of amplitude momenta? 

It was suggested in [15] to calculate not only the second, (jV'P)) and other 
momenta, but also physical values, such as (ln|?/'|). The average quantities 
are divided in [15] into two categories: containing physical information (e.g. 
(ln|V^|)) and containing no physical information (e.g. 

In other words, one has to understand how the choice of a semi-infinite system 
with a selected direction (n axis) and the boundary (layers in the transversal 
directions with n 0, 1) allows us to detect the localized states under question. 

Let us consider for simplicity a 1-D system, where the amplitude ip — ipn and 
the second initial condition i/^i = a. Again, for simplicity let us consider a to 
be a real quantity, i.e. [-^I^ = -0^. The recursive equation reads: 



where {Sn) = 0, {el) = a^. 

Let us consider now the random walk problem (Brownian motion [1,19] for 
the coordinate ipn) which is described by the following equations 



There exists a mean squared displacement (e^) = of the Brownian particle 



Ipn+l = (-B - en)lp„ - 1p„-l, 



(4) 



(5) 
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during a step. 



Considering n as the discrete-time index, both eqs. (4) and (5) describe the 
dynamics of the system with a stochastic time-dependent perturbation £„. 
The only difference is that the dynamics of the system (5) is triviah ip^ = ipQ 
when there is no perturbation, = 0. In contrast, the dynamical system (4) 
even for = reveals proper dynamics (it is a second order equation!). In 
the band region, \E\ < 2, this corresponds to the bounded motion. Therefore, 
the proper dynamics of both systems corresponds to the hounded trajectories. 



Assuming now that cr 7^ 0, eq. (5) describes the normal diffusion. The diffusion 
motion is characterized by divergences., e.g. for the mean time when the system 
returns to the initial state. In particular, the momenta of the amplitude 
are also divergent with the discrete-time n, e.g. 



reveals a power law divergence. For normal diffusion the mean square displace- 
ment ('i/^^) of the Brownian particles is linear in time. 

To detect the diffusion, it is sufficient to demonstrate the divergence of the 
second moment of the amplitude and to establish its law of time- dependence 
[20,21]. For example, for anomalous diffusion the mean quare displacement 
(t/'^) of the particles is no longer linear in time, but follows a generalized 
Fick's second law [22]. 

Generally speaking, other moments contain additional information, which, 
however, is not important. The divergence of the second moment defines the 
conditions of the diffusion appearance: if the second moment is divergent, so 
are all other even moments (this is why the choice of a particular moment 
for further analysis is not unique). The law of the divergence of the time- 
dependence allows us to distinguish between normal and abnormal diffusion. 
Important is the fact that the second moment (V^^) in eq. (5) can be found 
exactly analytically and thus we prefer its use (■0^-definition) . Speaking for- 
mally, the diffusion could also be classified using other averaged quantities, 
e.g. (|'0n|)) but such a choice is not convenient for mathematical reasons; it 
hinders an analytical solution. 

As is well-known, in a 1-D system described by eq. (4) any disorder cr 7^ leads 
to localization. This manifests itself in the simultaneous divergence, as n — > 00, 

of different average quantities of \iljn\^ e.g. the linear divergence of (ln|?/'„|) 
(log-definition of localization), whereas the exponential divergence occurs for 
the powers of \ipn\ [17]. The appearance of the localization in the approach 
based on eq. (4) is equivalent to the appearance of diffusion. In eq. (5) the 
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Fig. 1. Exponential divergence of second moment for generalized diffusion, eq. (4), 
for a = 0.1 and E = 0. 

random perturbation e„ is additive, which determines the linear (power-law) 
character of the second moment divergence (normal diffusion). In contrast, eq. 
(4) contains the product of and ipn {multiplicative noise) which determines 
the exponential character of the divergence (see also Brownian motion [3] and 
normal nonlinear diffusion [23] with multiplicative noise). 

I.e., we can speak of the generalized diffusion, analyzing the {ip'^) — f{n) 
divergence as a function of n. The exponential localization here corresponds 
to the exponential growth of f{n) oc exp(27n). Fig. 1, where ^ = I/7 can be 
interpreted as the localization length {ip"^ - localization definition [12,13,17]). 
On the one hand, the numerical solution for the amplitude ipn are plotted 
here using eq. (4) together with arbitrary realizations of perturbations Sn (two 
upper windows). On the other hand, the lower window of Fig.l show the exact 
analytical solution for the second momenta (ip^) (eqs.(51) and (53) of this 
paper). The analysis of these results (as well as other data for different E and 
a) demonstrates that purely exponential behaviour, {ip"^) oc exp(27n), is valid 
only asymptotically, as n — >^ cxd. In other words, the exponential behaviour 
of divergence is understood by its definition, eq.(49). There is a transient 
regime for finite n values especially for small 7 (as shown in Fig.l) where 
the exponential behaviour is combined with the oscillations corresponding to 
the proper dynamics, eq.(5). Note that the bounded proper motion of ip^ for 
£„ = in the band region, < 2, corresponds to oscillating functions. 

The ■0^-definition is convenient not only because the equations for the second 
moment can be exactly solved (see below). It is important that this defini- 
tion works also for the non-exponential localization, which corresponds to the 
non-exponential behaviour of /(n). In contrast, the traditional log-definition 
is valid only for the exponential localization, and does not allow an exact 
analytical solution. 
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It was analytically shown [17] for 1-D systems that if the higher moments of 
the random potential e„ can be neglected (for a 0) and we can restrict 
ourselves to the only parameter a, the exponents in all even moments (as 
well as in log-definition) are proportional to each other, they differ only by 
numerical factors. The proportionality of the different exponents (Lyapunov 
exponents) in the exponential growth indicates that different definitions of the 
generalized diffusion are in fact identical and one can use any of them, e.g the 
V'^-definition. 

Note that the emergence of the generalized diffusion phenomenon and the ex- 
istence of the phase diagrams in dynamical systems with the proper dynamics 
are analogous to critical phenomena in equilibrium and non-equilibrium sys- 
tems. In this sense, it is clear that use of the ■j/^^-definition for the determina- 
tion of the phase-diagram is quite sufficient. Exactly solvable problems in the 
physics of phase transitions [11], e.g. the Ising model, have demonstrated that 
in spite of the fact that the physical quantities, i.e. long range order parame- 
ter, magnetic susceptibility, specific heat, etc. are defined as different averages, 
they still behave similarly, when only the phase-diagram is considered. This 
is not surprising since one deals here with averages obtained for the same 
statistical ensemble. In particular, the critical temperature obtained for the 
order parameter coincides with that obtained for the susceptibility. This oc- 
curs for the exact solutions, although deviations can arise for the approximate 
solutions. Similarly, the phase diagram for the system with a metal-insulator 
transition can be defined uniquely if one determines any non-trivial average 
quantity. This is why the choice of (In [■^D is not convenient. A product of a 
sum of variables can be represented as a polynomial, where each term can be 
calculated, whereas the logarithm of a sum cannot be represented this way and 
thus the causality principle cannot be used. On the other hand, the averages 
of the powers of the amplitude can be calculated exactly. 

When we speak of the phase-diagram, we mean determination of the bound- 
aries of metallic and insulating phases. In this respect, the results [12,13] 
demonstrate that the second order moments form a closed and linear set of 
equations and, where the only characteristic of the random potentials is the 
second momentum, a^, have a simple interpretation. The phase boundaries 
are defined by the only simple parameter a of the random potential. That is, 
calculation of second moments (correlators) is the direct way to obtaining the 
phase diagram. 



2.3 Correlators 

In general, the field am is complex. We can introduce the simplest non-trivial 
correlators based on the ■0^-definition: 
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X{n)^^l = (CmV'n,l) , (7) 
yinU,l = ^[(C™^n-l,l) + (<-l,™^n,l)]- (8) 

Let us for mathematical simplicity consider real am- This is also justified by 
the fact that the final eqs. (36), (37) remain the same. For real ccm definitions 
can be simplified: 



x{n)mA = (^n,m^n,l) , (9) 
y{n)m,\ = {i>n,mi>n-l,l) ■ (10) 

Taking square of the both sides of the eq. (2) and using the causality principle, 
one gets 

x{n + l)m,i = Sm,icr'^x{n)m,in + x{n - + (H) 

^ra,ra'X{n)m',l'^l',l " £^m,m'y{n)m' ,1 " A,l'Z/(?T')l',m- 

Analogously one obtains the additional relation 

y{n + l)i„,i = -y{n\m + Cm,m'x{n)m',h (12) 

Note that cq.(ll) contains explicitly the diagonal correlators with 1 = m: 
x(n)in = x{n)m,m.- It is essential that diagonal correlators are positive num- 
bers, x(n)m > 0. 

2.4 Z-transform 



It is convenient to use along the n axis the Z-transform which is common in 
discrete-time systems [24] : 



= E ™> (13) 

n=0 

n^u = t (14) 

n=0 ^ 

After transformation one gets 



{Z- Z ^)X{z)m,\ - Cr'^Sm,\Xiz)m - X{l)m,\ = 
^m,m'-'^(^)m',l''Cl',l " ■^m,m'^(^)m',l — ^l^'Y ^m, 
zY{z)m,l = -Y{z\m + >Cm,m'^(^)m',l- 



(15) 

(16) 
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Hereafter the symmetry properties of the operator jCm,i are used. 
The solution of the second equation is simple: 



z 1 

Y(z)in,l = -^-—rj0.m,m'X(z)m',\ 2 _ 1 -^(^)m,l'>^l',l, (17) 



which leads to the equation for X-correlators 



[Z- Z ^)X{z)^^i - (T^5m,lX(2;)m " a;(l)m,l = j-^m,m'^(^)m',l' A',(ll8) 

Z j- 

~T [^m,m'^m', 

l'X(z)i>^l + ^(2;)i„,m''^m',l'A',l]- 



Z^ 



Taking into account the second initial condition, a;(l)in,i = amCKi- 
2.5 Fourier transform 



The obtained equation can be easily solved formally via Fourier expansion. 
Indeed, one gets for diagonal correlators 



xizU- /J^x(^,k)e-^'^™, (19) 
X(^,k)=5:x(^)me*". (20) 

m 

A general equation for X-correlators has to be solved using double Fourier 
expansion: 



X{z,Kk') = E^(^)m,ie^'^'"+''^'. (21) 

m,l 

This leads to 

U{z, k, k')X{z, k, k') = Q;(k)Q;(k') + a\{z, k + k'), (22) 
where 



C/(^,k,k') ^{z- z-^) - ^£(k)£(kO + -^[£^(k) +£2(k')], (23) 
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and 



a(k) = ^a^e^'^'". (24) 

m 



The function 

p 

>C(k) ^E-2Y, cos{kj) (25) 
arises as the Fourier transform of the operator >Cin,i- 

2.6 Diagonal correlators 

Assuming k" = k + k', eq.(22) can be rewritten as 

U{z, k, k" - k)X{z, k, k" - k) = a;(k)Q;(k" - k) + a\{z, k"). (26) 
Using the relation 

X{z,k") = j j^^X{z,k,k"-k), (27) 
eq.(26) can be solved elementary for diagonal correlators: 

Xiz,k)=niz,k)x^'\z,k), (28) 
v^^.k)- / ^'^^ c.(kOa(k-kO 

^ ^^'^^-y (2n)PU(z,k'k-k'y ^^^^ 



= 1 



(27r)P [/(^,k',k-k')^ 
2 r dPk' 1 

7^(z,k) ~ ^ J (27r)f C/(z, k', k - k') ■ ^ ^ 

Therefore we have obtained for an arbitrary second initial condition (field 

ctm) an exact solution for fundamental averages - diagonal correlators. After 
use of the Fourier transform the solution is expanded in modes labelled by 
the k-vector and describing the oscillations in the solutions in the transversal 
direction. It is easy to see that the V'^-definition leads to equations contain- 
ing only the a parameter of the random potentials. If cr = 0, all functions 
Tiiz.k) = 1 and thus xi^jk) = x'''^Kz,k). In other words, functions x^'^K^jk) 
correspond to solutions in a completely ordered system. Note that the bound- 
ary conditions (field a^) influence only functions x'^°)(2;, k) which are inde- 
pendent of the parameter a. Introduction of disorder, a > 0, transforms the 
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initial solution x^°H-2,k) into x(2;,k), where the operator 7^(2;, k) describes 
this transformation. 

In our derivation we assumed the system size in a transversal direction to 
be infinite, rrij e (—00,00), and thus the equations obtained correspond to 
the thermodynamic limit {L = 00). In the case of finite size of the system, 
rrij = 0, 1, . . . , L — 1 (cyclic conditions) the integrals should be replaced by 
series. Let us consider for illustration the 2-D case: p = D — 1 = 1. Here the 
index m. — m has L values and the number of modes k — k also equals L. 
That is, analysis of the diagonal correlators leads to L-order matrices (vec- 
tors). Surprisingly, this obvious fact was interpreted in ref. [15] that degrees 
of freedom of the initial problem were lost in our paper [12], thus question- 
ing our proof of the Anderson theorem. This conclusion was based on the 
traditional approach of the transfer- matrix [15,18] which operates with much 
larger, x L^, matrices for 2-D. The transition from x matrices to L- 
order matrices was interpreted as an error arising due to averaging over initial 
conditions [12] and thus artificially introducing a translation in the transversal 
directions. However, we do not use in the present study such an averaging (see 
also Section 3.6). The L-dimensional matrices are the natural mathematical 
tool for the description of diagonal correlators. Moreover, an excessive size of 
X matrices indicates that the transfer-matrix formalism is not adequate. 
This is seen, in particular, from the existence of the so-called trivial eigen- 
values of the transfer-matrix [15,18], which are cr- independent. Moreover, the 
transfer-matrix does not permit to make the transition to the thermodynamic 
limit and thus to calculate the phase-diagram (for more details see [16]). These 
disadvantages are absent in the 1-D case, when formally L = and the matrix 
dimensions coincide. This is why we used in [12] the transfer-matrix method 
only in the 1-D case and exclusively for the illustration. 



3 Anderson localization and stability 

3. 1 Signals and filters 

The set of equations for the correlators is linear, the same is true for eqs.(28) 
for each k-mode. The modes are normal (no mixing of modes with different 
k values). The above mentioned divergence (exponential or non-exponential 
growth) of the diagonal correlators, arising due to localized states, transforms 
mathematically to the instability of the set of linear equations. For linear 
discrete-time systems under study the most adequate formalism is signal the- 
ory [24]. Following [12,13,24], let us define x^^^^i'^) ^-s input signals. In our 
particular case this is a mathematical characteristic of the k-mode for an ini- 
tial ordered systems. (Its sense is explained below). The inverse Z-transform 
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of a given mode can be performed as follows: 



(31) 



The input signal is a one dimensional numerical series. The x{z, k) is the output 
signal; respectively after the inverse Z-transform the latter corresponds to the 
numerical series x(n, k), associated with the disordered system, a > 0. The 
relation between these two signals is given by cq.(28): the output signals are 
linear transformations of input signals performed through 7i{z, k) functions 
called the system filter [24]. 

Physical systems with time as independent variable are causal systems [24]. In 
our problem with discrete-time index n the causality is of primary importance 
for the interpretation of the result. Note that in the general case the inverse 
Z-transform is defined through the complex integral 



The integration here is performed over the complex plane called the region of 
convergence (ROC) [12,24], well-defined for the causal filters. For the causal 
systems h{n, k) = for n < holds, thus after the inverse Z-transform eqs.(28) 
transform into the convolution property [24] : 



One can see here the linear transformation of input signals into output signals. 

The fundamental property of the signal theory is that the divergence of the 
output signals is related to the asymptotic behaviour of the filter coefficients 
h{n, k) as n — > oo, which is mathematically equivalent to the poles of the 
7i(^,k) function of the complex argument z [24]. 

3.2 Fundamental mode 

Since the appearance of localized states leads to the divergence of the diagonal 
correlator, it is necessary to clarify, which fundamental mode k = ko looses 
its stability and thus is responsible for the divergence. Such a problem is quite 
general in many branches of physics. However, in our case it is quite obvious 
that the fundamental mode is static, ko = 0. 

The mode with k 7^ describes transversal oscillations of the diagonal cor- 
relators, but the sign of the divergent oscillating solution is not defined (see 




(32) 



n 




(33) 



1=0 
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also section 3.6). On the other hand, signs of the diagonal correlators are well 
defined, they are non-negative, x(n)m > 0. That is, in our particular case the 
solution is trivial, the fundamental mode is static. It is thus necessary to study 
the boundaries of the stability of the fundamental mode ko = only (the sign 
of the divergent non-oscillating solution is defined) , which means uniquely the 
determination of the phase diagram of the system. 

Assuming x^°Kz, 0) = S^°\z), x(^, 0) = S{z), H{z, 0) = H{z)), we arrive at 



S{z) = H{z)S^''\z). (34) 
After the inverse Z-transform one gets 



Sr. 



n = ^K-isf^- (35) 



The fundamental input signal S^^\z) here and the fundamental filter are de- 
fined by the integrals: 

S^^Hz) - / '^"^ l^^'^^l' (36) 

~ {z-l)J {2n)P [{z + \flz - £2(k)] ' ^""^^ 

1 . (^ + 1) / 1 

E{z) (^-l)i (27r)f [(^+l)V-2->C2(k)]" ^ ' 

Note also the relation 

h,,^^.lE{z)z^-. (38) 



27ri J z 

That is, we have got here a new derivation of relations earlier derived in Refs. 
[12,13]. Note that for the 1-D case the integral disappears, since p = 0, and 
one gets 



(39) 

1 -l-.^ii±^ I (40) 



H{z) {z-\)\{z^\flz-E-')\ 
3.3 Fundamental input signal 

Let us discuss the physical sense of the input signal signal S^^^z) (or s^"-*). 
Note that this corresponds to an ideal system with a — 0. When there is 
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no disorder, Em = 0, the particular solutions of the tight-binding cq. (1) 
are bounded functions, plane waves exp(i/C7W) with /C = {ki, k2, ■ ■ ■ , ko}, 
provided for a given energy E 

£{1C) = E (41) 



holds, where 

D 

£(/C) =^2cos(A;j), (42) 



otherwise particular solutions are not bounded functions, which lie beyond the 
band and have no physical interpretation. If we divide the wave vector into 
transversal or normal directions, /C = {k, ko}, and keeping in mind that 

C{k) = 2cos{kD), (43) 



for a given transversal mode k the bounded physical solution exists, provided 
|£(k)| < 2. These relations are sufficient for understanding the physical sense 
of the input signal. 



Eq. (36) contains in the integrand the function 

(^ + 1) 1 

(^-1) {z + iy/z- c^{k)' 

Under the condition |>C(k)| < 2 its inverse Z-transform gives 



(44) 



sin\kn) ' ^ ' 

where ko — koik.) is the solution of eq.(43). Assuming that the arbitrary field 
ccm is selected in such a way that the Fourier transform coefficients Q;(k) are 
not zero only if |>C(k)| < 2, eq.(36) after the Z-transform reads 



„(o) _ f |a(k)psin'(fcz)n) 

" J (27r)f sm^ko) ' ^ ^ 

The s^^ is bounded for any n. 

On the other hand, assuming that |>C(k)| < 2 is violated for certain k, so 
the a{k) ^ for |>C(k)| > 2, the inverse Z-transform of eq.(44) leads to the 
function 
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which increases without bounds as n 
of equation 

\£{k) \ =2cosh{KD). 

The input signal s^^-* is also increasing to infinity as a function of n. Note that 
both input and output signals are real values. 

In other words, on can conclude that any physical solution of an ideal sys- 
tem {<T = 0) with bounded wave functions has a one-to-one correspondence 
to the mathematical object - a bounded 1-D input signal Such solutions 
always exist, if the energy E lies within the band interval. On the other hand, 
formal solutions without physical interpretation {unbounded solutions beyond 
the band) correspond to the 1-D unbounded input signals. Note that partic- 
ular numerical values of input or output signals are not important from the 
point of view of the signal theory [24] whose main concern is a qualitative 
discrimination between bounded and unbounded signals. 

3.4 Fundamental filter 

The quantity characterizing the phase diagram of a disordered system is the 
fundamental filter, eq.(37). Its idea is quite clear [12,13]. Let us consider the 
physical states inside the band, \E\ < 2D. All these states are described 
by wave functions bounded in amplitude, which corresponds to a bounded 
1-D input signal s!^\ In a disordered system this signal transforms into s„, 
according to eq. (35). This output signal s„ can be either bounded, or infinitely 
growing (generalized diffusion). 

The natural interpretation of bounded output signals is that they - as before 
- correspond to the delocalized states bounded in amplitude. The unbounded 
output signals correspond to the localized states, respectively, which lead in 
the semi-infinite system to the divergence of the diagonal correlators. An im- 
portant result of signal theory [24] is that the divergence does not result from 
properties of a particular bounded input signal; the cause lies in the filter, 
H{z) or hn- We define the phase diagram of the Anderson model based on a 
general concept of the signal theory known as a BIBO stability [24]. Namely, 
a system is BIBO stable if every Bounded Input leads to a Bounded Output. 
A stable system (delocalized states) is characterized by a stable filter H[z); 
its main property is the absence of poles in the complex 2;-plane outside the 
circle \z\ > 1. An unstable system (localized states) is characterized by the 



(47) 

00. Here kd — /^^(k) is the solution 

(48) 
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unstable filter H{z) with poles outside the circle \z\ > 1. As was shown in 
[12,13], in the Anderson model the pole Zq of the unstable filter lies on the 
real axis. Consequently, Zq = exp(27), where 7 is the Lyapunov exponent 
(^^-definition) [12,13]. Using the inverse Z-transform, one can easily find that 
asymptotically /i„ ~ exp(27n), i.e. it is divergent. The localization length is 
defined as ^ = 7~-'^. 

The mathematical formalism of the signal theory allows to extend and comple- 
ment this result for the energy range \E\ > 2D, where in the absence of disor- 
der there existed only mathematical (divergent) solutions having no physical 
interpretation. It is known that disorder extends the band, i.e. new physical 
states arise also at \E\ > 2D. The question is: can delocalized states exist 
amongst these new states? The answer is simple. In the region \E\ > 2D only 
unbounded input signals exist. Such a signal cannot be transformed into 
a bounded output signal s„. The output signal s„ is always unbounded with a 
dual interpretation: on the one hand it corresponds to mathematical solutions 
(no physical interpretation); on the other hand, the divergence of the output 
signal can be associated with an emergence of the localized states. I.e., an 
emergence of delocalized states in the region \E\ > 2D (outside the band) is 
impossible [13]. 

Use of the filter function H{z) is a general and abstract method for describing 
the metal-insulator transition, valid for any space dimension. Earlier [12] we 
used a less fundamental approach for the 1-D problem, based on the calculation 
of the Lyapunov exponent for the ■0^-definition 

7= lim — lns„, (49) 



where for the 1-D case = (ipn)- Taking into account that the divergence of 
the output signal s„ is caused by the divergence of the filter one can use 
also the following equation 

7 = lim ^ In hn- (50) 



In the 1-D case both definitions of 7 are equivalent, since for a fixed energy 
and disorder, for a second initial condition ipi — a, one gets the single output 
signal Sn- 

Let us use for illustration eqs.(39), (40). Restricting ourselves by the band 
center E = (which simplifies the equations) and making the inverse Z- 
transform, one gets 
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"n 












.n 



(-lAo)l, 



(-lAo)l, 



(52) 



(51) 



2 



(53) 



provided = |ap[l — (— 1)"]/2 (bounded input signal). Indeed, simultaneous 
divergence of the output signal Sn and the filter hn arises due to the fact 
that the asymptotic behaviour of both quantities is determined by the same 
parameter, zq — exp(27) > 1. 

However, this equivalence is no longer valid for a space dimension higher than 
one. Since the second initial condition is defined by the field this cor- 
responds to a continuum of input signals s^lf> and, respectively, a continuum 
of output signals s„. That is, the definition (49) is not valid, since it is not 
clear that different signals should correspond to the same Lyapunov exponent 
7. Formally, a whole continuum of solutions should be analyzed. However, 
eq.(35) demonstrates that the fundamental Lyapunov exponent 7 does not 
depend at all on the field am, it is sufficient to define the Lyapunov exponent 
using only eq.(50). 

3.5 Phase diagram and multiplicity of solutions 

A more detailed analysis [12,13,16] reveals another problem of the definition 
(49): this is valid only in the 1-D case where for any disorder a only the phase 
of the locahzed solutions exists. 

As it was mentioned above, for the calculation of the fundamental filter /i„ 
using the inverse Z-transform, eq.(38), the contour integration over the so- 
called region of convergence (ROC) [24] is necessary, provided the filter under 
consideration is causal. It is shown [12,13] that dependent on the energy E 
and disorder cr, two general cases are possible for the Anderson model with 
D > 2. In the first case, the ROC hes outside the circle \z\ > 1. The filter 
hn is thus unstable and describes the localized states (insulating phase). In 
the second case, the ROC consists of two domains in the complex plane: one 
domain is inside the circle \z\ < 1, another one - outside this circle, \z\ > 1. 
Respectively, there are two ways to calculate the integral using eq. (38). The 
double solution arises, one describing delocalized states (filter hl^\ metallic 
phase) , the other locahzed states (filter h'^^ , insulating phase) . 

Consequently, in this range of parameters E and a the two phases can co-exist, 
this is why the metal-insulator transition in the Anderson model has to be an- 
alyzed in terms of first-order phase transition theory [12,13,16]. It is known 
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that complex (first or higher order) phase transitions introduced purely by 
noise can be generated for spatially extended systems [4]. For arbitrary ran- 
dom potentials the wave function can be either localized or delocalized (no 
co-existence!). However, in terms of statistics of an ensemble of random poten- 
tials, both localized and delocalized solutions can arise with a comparable prob- 
ability. Namely this comparability of probabilities is the main characteristic 
of the first-order phase transition, in contrast to the second-order transitions 
where either a pure metallic phase (no localized states) or a pure insulating 
phase (no delocalized states) should exist: co-existence is impossible. 



3 




Fig. 2. Phase diagram for 2-D case. 

For the illustration let us discuss here the phase diagram in terms of energy 
E and disorder a for 2-D case as shown in Fig. 2 (the band region \E\ < 4). 
The other dimensions are discussed in Ref. [13]. The diagram in Fig. 2 shows 
that all states with energies l-E"! > 2 are localized at arbitrarily weak disorder. 
The localization at the finite disorder parameter a > cro(-E), where according 
to Ref. [12] 

ao{E)=2{l-E'/4Y/\ (54) 



is also observed in the energy range \E\ < 2. In other words, only pure insu- 
lating phase (no delocalized states) takes place in these cases. The localized 
states are observed also in the rest of the phase diagram, namely, for small en- 
ergies, \E\ < 2, and small disorder, a < ctq^E). However, as it was mentioned 
above, the localized states could be observed only for a certain realization of 
the random potentials. For other realizations the delocalized states occur in 
the same range of energy and disorder, < 2 and a < ao{E). 

Phase co-existence creates a serious problem of the choice of an adequate math- 
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ematical formalism. Strictly speaking, when calculating the average quanti- 
ties on the ensemble of different realizations of random potentials, the contri- 
butions from pure metallic and insulating phases are considered as equivalent. 
However, such heterophase averages have no physical sense [12,13,16]. In a 
two-phase system with first-order phase transition one is interested in prop- 
erties of pure phases. From this point of view, the earlier introduced output 
signals s„ are also heterophase averages. This is why the derivation of equa- 
tions for signals is not a goal, but the tool for obtaining a more fundamental 
property - the filter H{z) (or As was noted [12,13,16], it is such a filter 
which reveals the multiplicity of solutions and which permits us to determine 
the phase diagram for the Anderson model. 

3.6 Averaging over initial conditions 

Let us discuss now the problem of averaging over initial conditions [12,13]. 
This procedure was introduced in Ref. [12], its meaning is quite simple. First 
of all, we choose a particular second initial condition, i.e. a field ccm- Secondly, 
let us consider the field a'^ = ttm+mo obtained from the first field as the 
result of a trivial translation in transversal direction by the vector mo. It is 
obvious that the relevant diagonal correlators x'(n)in = x(^^)m+mo can also be 
obtained by the argument shift in a transversal direction. Both solutions are 
physically equivalent. The Fourier transform of the field with the vector 
shift mo satisfies a simple relation: Q;'(k) = exp(— ikmo)Q;(k). 

Eqs. (28) to (30) clearly demonstrate that the vector shift afi'ects only signals 
for particular modes k, and all physically equivalent solutions differ from each 
other only by phases: 



Note that only the fundamental mode k = remains invariant. Averaging now 
the signal over all translations in transversal directions (averaging over initial 
conditions [12]) gives zero for all non-fundamental modes, k 7^ 0: x*^°^(-2, k) — 
0, x(2;, k) = 0. Nonzero is only the fundamental mode (k = 0): x(°)(2;, 0) = 
x(-2, 0) = S{z). The average diagonal correlator loses its dependence 
on the argument m: x(n)m = That is, averaging over initial conditions 
is an efficient tool for getting rid of all non-fundamental modes which are 
non-essential for the analysis. 

Note that eq. (36) contains r(k) = |Q;(k)p which is a function of the field c^m, 
its Fourier transform is T^. The averaging over initial conditions corresponds 
mathematically in the initial eq. (18) to the replacement of the initial condition 



X 



:'W(z,k) = exp(-ikmo)x(°n^,k), 
x'iz, k) = exp(-ikmo)x(^, k). 



(55) 
(56) 
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x{l)m,i = amtti by the average, x(l)in,i = ^m-i- After such a replacem^ent of 
initial conditions the system becomes much simpler: it is translation-invariant 
in transversal directions which makes use of the double Foiuicr transform 
unnecessary. In other words, averaging over initial conditions is nothing but 
a simple mathematical trick which permits to get quickly the fundamental 
property of system - the filter H{z). Such tricks are based on the fact that the 
filters are universal system characteristics, describing the solution transfor- 
mation from order to disorder. Such universal characteristics are independent 
of the initial conditions and other details. The initial conditions determine 
such non-universal properties as signals. Consequently, in order to obtain the 
fundamental filter H{z) and then the phase diagram, one can perform linear 
operations with signals, in particular, averaging over initial conditions. 

Another example of such operations is the following. Let us consider the field 
ctm not as fixed (second initial condition), but as a random variable, with si- 
multaneous averaging over the field and over the ensemble of the random 
potential realizations Sm- We assume the absence of correlations between the 
potentials Sm and the field unlike the correlation between the field com- 
ponents characterized by the arbitrary correlation function (amCti) = Fin-i. 
The initial condition in eq.(18) becomes a:(l)„i,i = ^ra-\- Taking into account 
that the operator Cm,m' by its definition, eq.(3), depends only on the on the 
argument difference, m — m', one gets the translation invariant system, where 
X{z)xii,\ = ^{z)m-i, provided for the diagonal correlators = x(^)m = 

X{z)o = S{z) holds (the diagonal correlators are m- independent). 

In this case in order to solve eq. (18), it is sufficient to use a single Fourier 
transform 



X{zU = J -—X{z,k)e-'^"^. (57) 



Instead of eq.(22) one gets 



U{z, k, -\s.)X{z, k) = F(k) + a'^S{z), (58) 
which gives 



^(-) = /(^^(-'k). (59) 

As a result, one returns to the fundamental eqs.(34) to (37), with replacement 
of the |a(k)|2 for F(k). 



21 



4 Conclusion 



We would like to stress that in this paper we presented a mathematically 
rigorous method for the calculation of the phase diagram for the Anderson lo- 
calization in arbitrary dimensions, which was briefly discussed earlier [12,13]. 
The phase diagram for the metal-insulator transition is obtained using the 
Lyapunov exponent 7. Localized states correspond to values of 7 > 0, i.e. 
a divergence of the averages over wavef unctions. This divergence is mathe- 
matically similar to the divergence of averages for the diffusion motion. That 
is, transition to the localized states can be treated as a generalized diffusion. 
From this viewpoint, in order to determine the range of the existence of local- 
ized states (i.e. the phase diagram) and the type of localization (exponential 
or non-exponential), it is sufficient to solve equations for the joint correlators. 
We have shown that these equations are exactly solvable analytically. 

In its turn, the appearance of the generalized diffusion arises due to the in- 
stability of a fundamental mode corresponding to correlators. The generalized 
diffusion can be described in terms of signal theory, which operates with the 
concepts of input and output signals and the filter function. Dclocalizcd states 
correspond to bounded output signals, and localized states to unbounded out- 
put signals, respectively. Transition from bounded to unbounded signals is 
defined uniquely by the filter function H{z), or more precisely, by the position 
of its poles in the complex plane. This function can be calculated for arbitrary 
space dimension D. 
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